Quantum spin configurations in Tb 2 Ti 2 07 
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Low energy collective angular momentum states of the Tb 3+ ions in Tb2Ti207 are classified 
according to the irreducible representations of the octahedral point group. Degeneracy lifting due to 
the exchange interaction is discussed. Diffuse neutron scattering intensity patterns are calculated for 
each collective angular momentum state and the ground state is inferred by comparing to experiment. 
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Tb2Ti207 is geometrically frustrated due to its py- 
rochlore crystal structure, apparently leading to a 
"spin liquid" state which persists to extremely low 
temperatures.— The spin liquid is characterised by an 
abundance of low-lying excited states which tend to pre- 
vent the formation of an ordered state because of fluctu- 
ations. The focus of experimental work has been to de- 
termine the conditions under which the spin liquid state 
becomes unstable against the development of long range 
order. Theoretical studies of antiferromagnetic interac- 
tions predict that ordering should occur with a Neel tem- 
perature of approximately 1.2 K)2j2i4 but this is not seen 
experimentally. There are reports of short-range ordering 
occurring at T c = 0.37 K— & and a spin glass phase below 
200 mKjSiZiS however one group has found that the spin 
liquid phase persists down to at least 50 mK. 1 ' 9 Even in 
the absence of long range order, one may expect that local 
correlated spin configurations, possibly governed by the 
exchange interaction, will be present. Indeed, a diffuse 
neutron scattering experiment that found an intensity 
pattern characterised by a broad peak at q = [0,0, 2] 10 
has been interpreted in this way?^ 

There is a vast literature of theoretical works on the 
general topic of spin liquids in geometrically frustrated 
systems A small subset of these treat the magnetic ions 
as classical Ising spins which point toward the corners of 
the tetrahedra, and this approach has been useful for un- 
derstanding spin ice behaviour in the related compounds 
Ho 2 Ti 2 07 and Dy 2 Ti 2 07i&i2 In this article, we extend 
this approach using basic considerations of symmetry and 
quantum mechanics. 

Tb 2 Ti 2 07 crystallises in the pyrochlore structure with 
space group Fd3m (#227, 7 h ). With two copies of the 
chemical formula per unit cell, there are four Tb ions 
which occupy the corners of a corner-shared tetrahedral 
network. The local site symmetry is 3m (D^d), where the 
three-fold axes point along the directions [111], [—1 — 11], 
[—11 — 1] and [1 — 1 — 1] for sites #1, 2, 3 and 4, respec- 
tively. In this environment, the J = 6 manifold of the 41* 
state of Tb 3+ splits into levels labeled by the irreducible 
representations (IR's) of D 3 : 3Ai © 2A 2 © AE^- Neu- 
tron scattering measurements and ab initio calculations 
find that the ground state and first excited states are E 
doublets, separated by 1.5 meV. 2 

Both the ground state and first excited state are lin- 



ear combinations of the angular momentum states | ± 5) , 
| ± 2), | =F 1) and | =p 4), where the quantum numbers 
measure J in the direction of the C3 axis. The ground 
state doublet was determined to be^iISi 



|±) = ±0.13| ± 5) T 0.13| T 1) - -97| =p 4). 



(1) 



With four Tb ions per unit cell, there is therefore a 
sixteen-fold degeneracy of the ground state for each unit 
cell. The collective angular momentum states are de- 
noted I ± ± ± ±) ee |±)i <g> |±) 2 <8> |±) 3 <g> |±) 4 , where 
the subscripts indicate the site number. These states are 
divided according to the representations by which they 
transform under the operations of the octahedral point 
group Oh, the underlying point group symmetry of the 
crystal, as A lg © 3E g © 2T lg © T 2g ~ The actual states 
are listed in Table I. The subscript g will be henceforth 
omitted since we will always refer to even representations. 

We begin by considering the nearest neighbour antifer- 
romagnetic (AF) exchange interaction. J^2 <ab> J a ■ Jb 
is isotropic, therefore it reduces the degeneracy of the 16- 
fold collective angular momentum states without mixing 
the different representations. To calculate the matrix el- 
ements, J a must be expressed in terms of the angular 
momentum quantisation axes (the C3-axis of each Tb 
ion). In the following, J a ; will always refer to the i-axis 
(i = x, y, z) using the local coordinate axes of the a posi- 
tion (a — 1,2,3,4), while J a (used below, in the discus- 
sion of neutron scattering) refers to global axes. For the 
specific choice of axes for the 1 and 2 positions, 



si = (l,l,-2)/V6, £ 2 = (-W, 
ft = (-1,1,0)/V2, y 2 = (1,-1,0 



(l,l,l)/\/3, z 2 



2)/V6 
2 

(-l,-l,l)/\/3, (2) 



one finds 



V2 



[Jlz{Jl+ + J2-) + [Jl+ + Jl-)Jlz] 



(Jl+J 2+ + JI-J2-) 
(Jl+J 2 _ + J1-J2+). 
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TABLE I: Basis functions of the collective angular momentum states for the four Tb ion sites, labeled according to the 
irreducible representations of Oh by which they transform, e = exp(27ri/3). 



For the axes of the 3 and 4 positions, 

fi 3 = (-l,l,2)/V6, Xi = (1,-1,2)/V6 
y 3 = (l,l,0)/\/2, y 4 = (-1,-1,0)/^ 
z 3 = (-1,1,-1)/V3, z 4 = (1, -1, -1)/V3 (4) 

some non-trivial phases appear: 

- - 1 
■n • J 3 = —-^Jiz'hz 
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Jiz(eJ 3 + + eV 3 _) + (eJi+ + e 2 J x _) J 3z ] 



+i(e 2 Ji+J 3+ + eJi_J 3 _) 



(Ji+J 3 _ + Ji_J; 



3+; 



(5) 



J\ ■ Ji — — - Jlz J42 



— [J lz (e 2 J 4+ + eJ 4 _) + (e 2 Ji+ + eJi_)J 42 ] 



+ -(eJi+J 4+ +e 2 Ji_J 4 -) 



G 



(Ji + </ 4 _ + J1-J4 



(6) 



where e = exp(27ri/3). The expression for J^'Jz is similar 
to Ji • J 4 , etc. The matrix elements are 



(±|J,|±) = ±3 

(±|J±|T) = i 



(7) 

(8) 



For the states given in Eq.Q] j = — 3.69?i and t is exactly 
zero (this is why the spins are Ising-like). The E^- 1 ' states 
have the lowest energy —2Jj 2 , while the E^ 2 \ Tj and 



T2 states have zero energy, and the At, E^ and 
states have energy 2Jj 2 /?>. Note that any linear com- 
bination of Ai, E^ and resembles the "spin ice" 
state (two spins pointing in and two spins pointing out 
of each tetrahedron), which is the ground state for fer- 
romagnetic exchange on classical spins. The E± states 
resemble the classical states often described in the lit- 
erature as the spins pointing outward from the corners 
of the cube which contains the tetrahedron, which also 
minimise the AF exchange energy. 

Because t = 0, the exchange interaction leads to no 
mixing between states listed in Table I belonging to 
the same representation, although this would be per- 
mitted under general considerations of symmetry. Non- 
zero matrix elements that mix |+) and |— ) states could 
originate from mixing with higher crystal field levels^ 
or with higher multipole interactions, beginning with 
quadrupole. Obviously such effects are relevant if the ex- 
change interaction constant J is small. If t then the 
mixing between states due to the exchange interaction 
is easily found using the matrix elements given in Table 
III, below. The dipole-dipole interaction, when limited 
to a single tetrahedron, only renormalises the exchange 
interaction, with slight differences between the diagonal 
and off-diagonal matrix elements of the sixteen collective 
angular momentum states. 

To determine which of the sixteen states is the most 
likely ground state for Tb2Ti2C>7 we compare the diffuse 
neutron scattering intensity pattern with calculations of 
the scattering function. In the dipole approximation, the 
scattering intensity is proportional to^i 

7(q) oce- 2 ^)|F(q)| 2 ^ e -^/ fe ^/ m ( g ) (9) 
m 

where W(q) is the Debye- Waller factor, F(q) is the form 
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factor for Tb 3+ and 



/m(q)=EEE^-^)H4l»)<»|J6'He <q - (rb - r » ) 

i,j a.b n 

(10) 

a and b are the four sites at the corners of the tetrahe- 
dron and n and m are the sixteen states, whose energy 
differences are assumed to lie within the energy range 
across which the neutron scattering is integrated. The 
functions / m (q) were calculated for each state m and the 
results for states belonging to the same doublet or triplet 
were added. The cases j ^ and t ^ were considered 
separately. The calculations yield four different patterns, 
shown in Fig. 1. The correspondence between different 
patterns and states is summarised in Table UTl 



state 



Ai 

E w 

E (S) 

r (D 

r (2) 
T-2 



aj 
-a) 



b) 

b) 
-b) 
c) 

d) 



TABLE II: The first column lists the states given in Table I, 
and the second and third columns indicate to which pattern 
of Fig. 1 each state corresponds, for non-zero j (t = 0) and 
non-zero t (j = 0) respectively. A negative sign indicates 
that the intensity is reversed, and no entry indicates that the 
pattern has a uniform intensity. 

If the exchange interaction plays a dominant role in 
this system then one would expect that Fig. la), which 
is associated with the anti-ferromagnetic state E^ and 
the ferromagnetic states A±, E^ and for j ^ 0, 
would match diffuse neutron scattering measurements. 
However, this is clearly not the case. Neutron scattering 
patterns oscillate in three directions, [0, 0, Z] , [0,/,/] and 
[/,/,/], and show a broad peak at the position [0,0,2]^ 
The only figure which shows a qualitative agreement 
to experiment is Fig. lc), and in fact the agreement 
is very good. Considering that the experimental re- 
sults are a sum over all / m (q), weighted by the fac- 
tor exp(— .Em/feeT), we conclude that the configuration 
leading to Fig. lc) is the ground state. 

Fig. lc) comes from the state , which is neither 
ferro-magnetic nor anti-ferromagnetic and completely in- 
sensitive to the parameter j. This result is difficult to 
explain, but at least it is self-consistent. The interaction 
which could lead to a t[ x ^ type ground state is unknown, 
but obviously it must contain an effective non-zero t. 
Likewise, the neutron scattering pattern shown in Fig. 
lc) can only be derived with a non-zero t. Moreover, a 
non-zero t restores the true quantum mechanical nature 
of this system, for without it only classical Ising-like spins 

i-om^n 4.18 



The origin of non-zero t is still an open question but 
there are a couple of possibilities worth considering. 
First, higher order multipole interactions restore quan- 
tum effects since the matrix elements 



(-\JzJ-\+) =(-\J-J z \+) 



-5.217T 



and 



[-\Ji 



= -0.71a 2 



(ii) 



(12) 



are non-zero. Alternatively, extending the 16- 
dimensional manifold to include higher crystal electric 
field levels will add some non-zero matrix elements for 
the J± operators. This approach was followed in Ref. [ill 
by including the first excited crystal electric field doublet, 
thus increasing the number of collective angular momen- 
tum states from 16 to 256. In that case, the calculated 
diffuse neutron scattering intensity pattern agreed well 
with experiment but a singlet ground state was found 
based on exchange and dipole-dipole interactions. 

The effect of a non-zero t on the exchange interaction 
on the 16 states given in Table I can be calculated using 
© ([5]) and ©■ The results are summarised in Table [TTT1 
The lowest energy state is found in the E sector for any 
values of j and t. Therefore, introducing a non-zero t into 
the exchange interaction rearranges the energy levels but 
never produces a T\ ground state. 



Sector 


Exchange matrix 


t = 


3=0 


Ax 


2j 2 /3 - 2t 2 /3 


2j 2 /3 


-2t 2 /3 


E 


( -2j 2 t 2 \/2j3 \ 

-t 2 /2 -4tj/V3 
\t 2 i/2/3 -4tj/V3 2j 2 /3 + t 2 /3 J 


2f/3 


-2f 


t 2 
-t 2 /2 
-2t 2 /3 


Ti 


( -t 2 /2 \ 
V 2j 2 /3 ) 


2j 2 /3 





-t 2 /2 


T 2 


5t 2 /6 





5i 2 /6 



TABLE III: Matrix elements for the exchange interaction for 
each representation within the 16-state manifold spanned by 
the states in Table I. The last two columns give the eigenvalues 
of the matrices for t = and j = respectively. 

To summarise, we have classified the low energy collec- 
tive angular momentum states of the Tb 3+ ions according 
to the irreducible representations of the octahedral point 
group Oh- Diffuse neutron scattering patterns are best 
described by the ground state. This result cannot 
be accounted for by the exchange interaction. Then, two 
important questions remain, namely, what is the nature 
of the local interactions leading to such a ground state, 
and what are the low-lying excitations responsible for 
spin liquid behaviour. 
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FIG. 1: (Colour) Contour plots of diffuse neutron intensity patterns /m(q) fEg . 110(1 for different states for non-zero j (figure a) 
and non-zero t (figures b-d) (refer to Table II.) In all four figures the intensities have been scaled in order to draw them using 
the same colour range. 
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